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Measurement of 3 with Dalitz Plot Analysis of B + — ► Decay 

A. Garmash a 

a Department of Physics, Princeton University, Princeton, New Jersey 08544, USA 

Results from the Belle and BaBar experiments on measurement of the weak angle <j>3 using a Dalitz plot analysis 
of the Kgn + n~ decay of the neutral D meson from the B ± -> D M K M± process are presented. The method 
employs the interference between D° and D° to extract the angle 03, strong phase 5 and the ratio r of suppressed 
and allowed amplitudes. 



1. Introduction 

Determinations of the Cabibbo-Kobayashi- 
Maskawa matrix elements provide important 
checks on the consistency of the standard model 
and ways to search for new physics. The possibil- 
ity of observing direct CP violation in B — > DK 
decays was first discussed by I. Bigi, A. Carter 
and A. Sanda [1]. Since then, various methods 
to measure the weak angle 03 (also known as 
7) using B — > DK decays have been proposed. 
All these methods are based on two key observa- 
tions: neutral D° and D° mesons can decay to a 
common final state, and the decay B + — ► DK + 
can produce neutral D mesons of both flavors via 
b — > cus and b — ► ucs transitions, with a relative 
phase 9+ between the two interfering amplitudes 
that is the sum, <5+03, of strong and weak interac- 
tion phases. For the decay B~ — ► DK~, the rel- 
ative phase is = 5 — 03, so both phases can be 
extracted from measurements of such charge con- 
jugate B decay modes. (Unless stated otherwise 
charge conjugation is implied throughout this re- 
port.) However, the use of branching fractions 
alone requires additional information to obtain 
03. This is provided either by determining the 
branching fractions of decays to CP eigenstates 
(GLW method [2]) or by using different neutral 
D final states (ADS method [3]). 

A Dalitz plot analysis of a three-body final 
state of the D meson allows one to obtain all 
the information required for determination of 03 
in a single decay mode. Three body final states 
such as KgTr + ir~ have been suggested [4,5] as 
promising modes for the extraction of 03. In the 



Wolfcnstein parameterization of the CKM matrix 
elements, the weak parts of the amplitudes that 
contribute to the decay B + — > DK + are given 
AX 3 (for the D a K+ final state) and 
A\ 3 (p + irj) (for D°K+). The two am- 
plitudes interfere as the D° and D° mesons decay 
into the same final state Kgir + ir~; the admixed 
state is denoted as D±. Assuming no CP asym- 
metry in neutral D decays, the amplitude of the 
D± decay as a function of Dalitz plot variables 
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where /(mq_,TO_) is the amplitude of the 
D®^K^-k + ti~ decay, and r is the ratio of the 
magnitudes of the two interfering amplitudes. 
The value of r is given by the ratio of the CKM 
matrix elements |^KJ/IK*f>KJ an d the c °l° r 
suppression factor, and is estimated to be in the 
range 0.1-0.2 [6]. 

The method has a two- fold ambiguity: the 
(03, S) and (03 + 180°, 8 + 180°) solutions cannot 
be separated. The solution with < 03 < 180° 
is chosen. 

The method described above can be applied to 
other B decay modes such as B + — > D*K + and 
B+ -> DK*+. 



2. Data Samples 

Results from the two B- factories Belle/KEKB 
and BaBar/PEPII arc available. The current 
proceedings are based on results reported in 
Refs. [7,8]. For the most recent updates see 
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Table 1 



Fit results for D°- 


-*KgTr + Tr~ decay [7]. 


Charmel 


Phase (°) 


Fit fraction 


K^i 


212 ±3 


9.8% 


K° sP Q 


(fixed) 


21.6% 




110.8 ± 1.6 


0.4% 


Kg/o(980) 


201.9 ± 1.9 


4.9% 


K° s a 2 


237 ± 11 


0.6% 


Kg/ 2 (1270) 


348 ±2 


1.5% 


tfg/o(1370) 


82 ±6 


1.1% 


K%p°(U50) 


9±8 


0.4% 


K*(892)+n- 


132.1 ±0.5 


61.2% 


K*{892)-ir+ 


320.3 ± 1.5 


0.55% 


iT*(1410) + 7r- 


113 ±4 


0.05% 


K*(U10)-tt+ 


254 ±5 


0.14% 


K*(U30)+ir- 


353.6 ± 1.2 


7.4% 


^(1430)-7T+ 


88 ±4 


0.43% 


isr 2 *(1430)+7r- 


318.7 ± 1.9 


2.2% 


k{{ 1430) -7T+ 


265 ±6 


0.09% 


if*(1680)+7r- 


103 ± 12 


0.36% 


if*(1680)-7r+ 


118 ±11 


0.11% 


non-resonant 


164 ±5 


9.7% 



Ref . [9] . The Belle collaboration uses a data sam- 
ple that consists of 386 x 10 6 BB pairs. The 
decay chains B+ -> DK+, B+ -> D*K+ with 
D*->Dn° and B+ -> DK*+ with K*+^K%tt + 
are selected for the analysis. Analysis by the 
BaBar collaboration is based on 227 x 10 6 BB 
pairs. The reconstructed final states are B + — » 
DK+axid B+ -> D*K+wiih two D* channels: 
D*^Dtt° and D*^Dj. The neutral D meson 
is reconstructed in the K^it + it~ final state in all 
cases. The dominant backgrounds come from a 
random combination of a real or fake D^*'° me- 
son with a charged track in continuum events; 
from B + — > D^it + events with n — K misiden- 
tification or other BB decays. 

The numbers of reconstructed signal events in 
the Belle's sample are 331 ±23, 81 ± 11 and 54±8 
for the B+ -> DK+, B+ -> D*K+ and B+ -> 
DiC* + channels, respectively. BaBar finds 282 ± 
20, 90 ± 11 and 44 ± 8 signal events in the B + -» 
DK+, B+ -> D*[Dtt°}K+ andB+ -> D*[D 1 ]K+ 
respectively. 

Several high statistics control samples are re- 



Figure 1. (a) m 2 _, (b) m 2 + , (c) and 
(d) Dalitz plot distribution for D*~ — > D°ir~, 
D°-^Kgir + ir~ decays from the e + e~ — > cc con- 
tinuum process. The points with error bars show 
the data; the smooth curve is the fit result [7]. 



constructed in order to cross-check the CP fit re- 
sults. A sample of D*~ — > D°tt~ events produced 
via the e + e~ — » cc continuum process is selected. 
This decay mimics the B~ D°K~ decay with 
r = 0. The B+ — > Z)(* ) 7r + mode is used as a 
second control sample where r is expected to be 
approximately 0.01. 



3. D° — >Kg7v + 7v decay amplitude 

The D°^Kgn + n~ decay amplitude 
/(m3_,m_) can be determined independently 
from a large sample of flavor-tagged D*~ — > 
D°it~ , Z)^ — > -^-s 7r^ 7r — decays produced in contin- 
uum e + e~ annihilation. Once / is known, a fit 
to B^ data allows determination of r, fa and S. 

The amplitude / is parametrized as a coherent 
sum of quasi-two-body amplitudes Aj and a non- 
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Figure 3. Results of signal fits with free pa- 
rameters x — r cos 9 and y — r sin for (a) 
B+ -> L>if +and (b) 5+ -> D*/v+ samples, sep- 
arately for f?~ and B + data. Contours indicate 
integer multiples of the standard deviation [7] . 



by minimizing the negative likelihood function 
(Belle performs an unbinned fit while BaBar uses 
a binned approach.) 



Figure 2. Dalitz plot and one dimensional pro- 
jections for the B~ — > D°K~ (left column) and 
B + — > D°K + candidates (right column) in the B 
signal region [8] . 



resonant amplitude, 
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E« 

3=1 
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where 2V is the total number of resonances and Oj , 
«j, 6 and /? are free parameters. A set of 18 quasi- 
two-body amplitudes is used to fit the data. The 
list of resonances included in the model as well as 
parameters determined from the fit are summa- 
rized in Table 1. For consistency with other an- 
alyzes [8,10], the D° — > Kgp° mode is chosen to 
have unit amplitude and zero relative phase. The 
D°^KgTr + TT~ Dalitz plot distribution, as well as 
its projections with the fit results superimposed, 
are shown in Fig. 1. 



4. CP Analysis 

Analysis of CP violation is performed by means 
of a maximum likelihood fit. The fit is performed 



-21ogL = -2^1ogp( 



2 2 



0,(3) 



with the Dalitz plot density p represented as 
p(m 2 + ,m 2 _,AE,M hc ) = 

e\f(m 2 + ,ml) + (x+iy)f(m 2 _,m 2 + )\ 2 F sig (AE,M hc ) 
+ F hck {m 2 +7 m 2 _,AE,M hc ) 7 (4) 



-^sig IS 



where x = rcos(03 +5), y — rsin(</> 3 +5) 
the signal AE and M^ c distribution represented 
by the product of two Gaussian functions (BaBar 
also includes Fisher discriminant in the fit); -F^ck 
is the distribution of the background events, and 
£ = e{m 2 L,rn?_) is reconstruction efficiency de- 
termined from MC simulation. The background 
density function Fbck is determined from analysis 
of sideband events in data and with Monte Carlo 
generated events. 

The fit procedure is first tested on two high 
statistics control samples: D* + ^D°tt + from cc 
continuum events and B~ — hDM 7r~ . In the CP 
fit to D*+^D Q 7r+ the fit yields r B = (-5.2 ± 
5.2) x 10~ 3 . Results consistent with r = 
are obtained with the other two control samples: 
r B = (1.8 ± 1.5) x 1(T 2 for the B-^D^tt- mode 
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Table 2 

Summary of fits results. 



Parameter 


Belle 


BaBar 


4>3 


53.3^7? ±2.5 ±8.7 


70 ± 311™ 


TDK 


0.159±g;g|^"=b 0.012 ± 0.049 


0.12 ±0.08 ±0.03 ±0.04 


5dk 


145.7+^;? ±3.0 ±22.9 


104 ± 45±%±% 


td*k 


0.175t£o99± 0.013 ± 0.049 


0.17 ±0.10 ±0.03 ±0.03 


Sd'K 


302.0+^ f ± 6.1 ± 22.9 


296 ±411^ ± 15 


TDK* 


0.564;°;^' ±0.041 ±0.084 




6dk* 


242.6+^ j ±2.5 ±49.3 





andr D . x = (4.6±2.1) x 10~ 2 for the B-^D*°tt- 
mode. This is in agreement with expectations. 

A summary of the results of the fits to the sig- 
nal events is given in Table 2, where the first 
quoted error is statistical, the second is system- 
atic and the third is a model uncertainty. Figure 3 
demonstrates results of the fit to B + — > DK + and 
B + — > D*K + events on the x — y plane. 

Systematic errors come from uncertainty in the 
knowledge of the functions used in the signal 
Dalitz plot fit. These include the Dalitz plot pro- 
files of the backgrounds and the detection effi- 
ciency, the momentum resolution description, and 
on the parametrization of the AE — Mt, c shape of 
the signal and background. Though the statisti- 
cal error is still quite large, this method currently 
provides the best direct constraints on $3. 

Uncertainty in the model used to parametrize 
the D°^Kgir + ir~ decay amplitude is the source 
of the associated error in the analysis. It 
arises from a non-unique choice of the set of 
quasi-two-body channels as well as uncertainty 
in parametrization of some components (non- 
resonant amplitude, for example). To evaluate 
this uncertainty several alternative models have 
been used to fit the data. 

Although at present the 03 accuracy is dom- 
inated by the statistical uncertainty, the model 
error will eventually dominate as the experimen- 
tal statistics increase. A model independent way 
to extract ^3 has been proposed in Rcf. [4]. 
The idea is that in addition to flavor tagged 
D°^KgTT + TT~ decays, one can use CP tagged 
decays to K^ir+ir- from the ijj(3770) DD 
process. Combining the two data sets the am- 



plitude and phase could be measured for each 
point on the Dalitz plot in a model independent 
way. Study with Monte Carlo simulations indi- 
cates that with with 50 ab _1 of data ^3 can be 
measured with a total accuracy of better than 2 
degrees. 
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